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ABSTRACT 



Context. The /3 Pictoris system with its debris disk and a massive giant planet orbiting at 9 AU represents an ideal laboratory to 
study giant planet formation and evolution as well as planet-disk interactions. yS Pic b can also help testing brightness-mass relations 
at young ages. Other planets, yet undetected, may of course be present in the system. 

Aims. We aim at putting direct constrains on the mass of yS Pic b and at searching for additional jovian planets on orbits closer than 
typically 2 AU. 

Methods. We use high precision Harps data collected over 8 years since 2003 to measure and analyse p Pic radial velocities. 
Results. We show that the true mass of /? Pic b is less than 10, 12, 15.5, 20 and 25 MJup if orbiting respectively at 8, 9, 10, 11 and 
12 AU. This is the first direct constraint on the mass of an imaged planet. The upper mass found is well in the range predicted by 
brightness-mass relations provided by current "hot start" models. We also exclude the presence of giant planets more massive than 
2.5 MJup with periods less than 100 days (hot Jupiters), more massive than 9 MJup for periods in the range 100-500 days. In the 
500-1000 day range, the detection limit is in the brown dwarf domain. Beyond the intrinsic interest for y6 Pic , these results show the 
possibilities of precise RV measurements of early type, rapidly rotating stars. 

Key words, stars: early-type - stars: planetary systems - stars: individual (J3 Pic ) 



1. Introduction 

Since the imaging of a circumstellar debris disk in the eight- 
ies, the youn g (12^^ Myr; IZuckerman et al. 2001) 1 and close 
(19.3^0.2 pc; ICrifo et al. 1997] l A5V star /3 Pictoris has been 
considered as a prototype of young planetary systems. Together 
with the other debris disks imaged since, this system allows to 
study the physical and chemical characteristics of sites of on- 
going or just finished planetary formation. Recently, we were 
able to detect a companion orbiting the star at a distance rang- 
ing between 8 and 15 AU (Lagr ange et al. 2009 ). Its U - 11.2 
approximate apparent magnitude translates into a temperature 
of ~ 1500 K and a mass of ~ 8 Mjup according to Lyon's 
group models dBaraffe et al. 2003T l. Similar masses were later 
found from images at Ks ( [Bonnefoy et al. 20111 , and at 4 
( [Quanz et al. 2010) 1. If orbiting on a slightly inclined orbit, the 
companion could explain most morphological (asymmetries) 
and dynamical peculiarities of the p Pictoris dust system, as 
well as the "falling evaporating bodies" (FEB) phenomenon at 
the origin of the replenishment of the disk gas phase. In par- 
ticular, it could explain the characteristics of the warp observed 
within the inner part of the disk (Lagrange et al. 2010 and ref. 
therein). Additional data now show that the planet semi-major 
axis is probably in the range 8-12 AU (see IChauvin et al. 201 II 
and references there-in). Interestingly in the present context 
of planet formation theories, jS Pic b is located in the re- 
gion where, according to [Kennedy and Kenyon 2008 models, 
and given the star age and mass, giant planets can have 
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formed in situ by core-accretion, in contrast to the few other 
planets imaged so far dMarois et al. 20081 IMarois e t al. 20101 
Kala^etal. 20081 ILafreniere et al. 20T0l IChauvin et al. 20051] 
iChauvin et al. 2005bl) . located further away from their parents. 

The mass determination of p Pic b as well as those of 
all other planets imaged relies however on the model depen- 
dant brightness-mass relation provided by the so-called "hot- 
start" models. In these models, the giant planet is formed by 
the collapse of a gaseous cloud and all the energy released is 
transfered in heat. The young planet is therefore quite hot dur- 
ing its first 100 Myr, which may allow planet direct detection 
around young stars with current imagers. Recently, however, 
[Fortney et al. 2008 have developped an alternative model, ex- 
pected to provide brightness-mass relationships in the case of 
core-accretion. In their model ("cold-start" model), a signifi- 
cant amount of energy is lost during the accretion process of 
gas onto a solid core. As a result of this loss of energy, their 
model predicts planets much fainter than the "hot-start" model 
predictions at young ages, at a given mass: for instance, a 10 
Jupiter mass planet aged 10 Myr is ^ 5-6 magnitudes fainter 
Under such assumptions, and except during the very short accre- 
tion phasis, young planets would not be detectable with current 
instruments. We note that the discrepancy between both mod- 
els decreases with time, but remains significant up to 100 Myr 
Coming back to the y6 Pic b case, the observed magnitude is com- 
patible with a planet formed under the "hot-start" model, but not 
under the "cold-start" model, as presently developped. In par- 
ticular, this "cold-start" model does not predict any planet with 
such luminosity at 10 Myr However, it would be prematurate to 
conclude that fi Pic b did not form by core-accretion based on 
these brightness-mass considerations, as these models, 1) prob- 
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ably represent two extreme situations in between which actual 
processes at work during giant planet formation take place, and, 
2) have not been calibrated so far by means of observations. 
Calibrations require independant measurements of the mass of 
the planets, ie dynamical masses, derived from astrometric or 
RV measurements. This is why detecting and characterizing 
planets, in particular young ones, is so badly needed. 

Another important issue is the presence of other planets 
around (i Pictoris. The example of HR8799 jMarois et al. 201 Oi l 
and the radial velocities (hereafter RV) study of solar-type 
stars (http://www.exoplanet.eu) show that multiple planetary 
systems are indeed very common. Also, gaps have been re- 
ported in the fi Pic disk in the mid-IR (Okamoto et al. 20041 
Wahhaj et al. 2003 ITelesco et al. 2005l l that were sometimes at- 
tributed to planets (IFreistetter et al. 20071 1. Our available images 
do not show evidence of planets more massive than ~ 5 M y„p at 
separations of 10-20 AU (Lagrange et al. 2010 1. As current im- 
agers are not able to put constrains on the presence of planets 
closer than typically 5-6 AU from the parent stars, other meth- 
ods are needed to get information on the presence of planets 
in the inner 5 AUs. An analysis of 120 spectra obtained with 
Coralie in 1998-1999 as well as 258 Harps spectra obtained be- 
tween November 2003 and July 2004 allowed us earlier to ex- 
clude the presence of an inner giant planet with masses larger 
than 2 My,,,, at a distance to the star of 0.05 AU, and 9 MJup at 
1 AU (|Galland et al. 2006). In this paper, we analyse more than 
1000 Harps spectra recorded between 2003 and 20 11 (described 
in Section 2) to search for additional, closer planets (Section 3) 
and to further constrain p Pic b mass (Section 4). 



2. Data 

We have now gathered 1098 Harps spectra from the ESO 
Archive or obtained during our observing runs. These data span 
a period of time from October 2003 to February 201 1 . Exposure 
times of individual spectra are between typically 60 and 120 sec- 
onds, and typical SNs vary between 150 at 250 at 400 nm. 

All data have been reduced and analysed with our SAFIR 
software dGalland et al. 2005) . that we developped to measure 
RV of early type stars precisely enough to detect planets. Such 
stars are usually rapid rotators and show much fewer lines than 
solar type stars; they are then not adapted to masks-based tech- 
nics. The method consists in correlating, in the Fourier space, 
each spectrum and a reference spectrum built by summing-up 
all the available spectra for this star to measure its RV. 

The RV thus obtained are shown in Figure [1] The most 
obvious feature is the presence of high frequency (HF) 
variations, already described in Coralie and Harps spectra 
dGalland et al. 2006 1) at 47.4^0.01 cycle/day (period of 30.4 min) 
and 39.05^0.01 cycle/day (period of 36.9 min) and part of a com- 
plex pattern of pulsations described in IKoen et al. 20031 Before 
2008, each pointing (refered to as a visit) consisted generally in 
two consecutive exposures, and lasted typically 5-10 minutes, 
including overheads. Starting from March 2008, we adopted an 
observational strategy aimed at attenuating these HF variations. 
The duration of the visits was increased to ^ 2 hours, during 
which we recorded several 1-2 minutes consecutive exposures, 
to average out as much as possible the HF variations. Examples 
of the RV obtained during these 2hour-visits are provided in 
Figure |2] The period March 2008 (JD 2454542) until February 
2011 (JD 2455597,referredtoas"setl"; 810 spectra over ^ 1000 
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Fig. 1. Top: RV data of jS Pic between 2003 and 2011. Bottom: 

associated periodogram. 
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Fig. 2. Examples of high frequency RV variations. 

Table 1. rms (m/s) of all RV measurements between 2003 and 
2011 and rms of the uncertainties (m/s) associated to the RV 
measurements. 



Data set 


no average 


average 0.5d 


average Id 


all data 


275.2 (39.1) 


257.7 (27.3) 


248.1 (27.3) 


setl 


277.6 (39.7) 


97.2 (5.2) 


97.2 (5.2) 



' Such masses measurements are refered in the following to as "di- 
rect", in contrast with model-dependent, brightness-based estimations. 



days) defines then a much more homogeneous set of data than 
the whole set of data. The impact of this strategy on the RV curve 
is illustrated in Figure [3] where we show the data once averaged 
over 0.5 day and 1 day. The improvement is also clear when we 
consider the rms of the data (see Table [U, which is significantly 
decreased once the data have been averaged. 

3. Search for short period planets 

3. 1 . rms-based detection limits 

Our SAFIR software provides estimations of the detection lim- 
its. These limits are computed in the following way: for a given 



2 



A.-M. Lagrange et al.: Constrains on planets around p Pic with Harps radial velocity data 



1,2 



0,8 - 




-2,8 - 



-3,2 -I 1 1 1 \ 1 1 

500 1000 1500 2000 2500 3000 3500 4000 

JD-2452000 

Fig. 3. RV data of fi Pic between 2003 and 201 1, averaged over 
0.5 (red) and 1 (green) day. For comparison, the unaveraged data 
are also given in blue. Averaged data have been vertically shifted 
for clarity purpose. 

planet (Mass, period), we compute the induced RV at the same 
observing dates as the real observations and the rms of these sim- 
ulated RV, for several (typically 1000) values of orbital phase. 
The distribution of the obtained rms is gaussian. The planet is 
detected if the rms of the observed RV (that we call hereafter 
RV-jitter) is smaller than the average of the rms distribution of 
the simulated RVs. The level of confidence (detection probabil- 
ity) is computed by comparing the standard deviation of the sim- 
ulated distribution and the gap between the observed rms and 
the average value of the simulated distribution. We showed in 
[Lagrange et al. 2009| that when the signal is well sampled, the 
detection limit is close to the one that corresponds, for a given 
period, to a planet that would produce an RV amplitude equal to 
3xRV-jitter This computation of the detection limits is very fast; 
however, it does not take into account the temporal structure of 
the stellar noise. As the fi Pic variations are mainly HF variations 
(Figure[TJ, this method tends to overestimate the detection limits 
for periods larger than one day (unless, when the data sampling 
allows it, temporal averaging can be done). We therefore devel- 
oped in addition alternative methods that take into account the 
frequency-structure of the stellar noise. These approaches will 
be described in Sect. 3.2. 

The obtained detection limits (99.7% probability) as com- 
puted on the whole set of data are provided in Figure H] Typical 
values obtained for P ^ 5, 10, 100, 500 and 1000 days are resp. 
7, 10, 16, 62 and 70 M y„p (all the computed masses given in this 
sub-section are rounded to the next integer or integer/2 mass). 
Note that as the system is seen edge-on, within at most a few de- 
grees, if we assume that the planets orbit within the debris disk, 
then the detection limits derived are not affected by the sin(i) 
indetermination (i=90 degrees). 

When averaging the RV over one day, we find detection lim- 
its of 4.5, 5.5, 13.5, 29 and 53 Mj„p (see also Figuregl. Hence, 
averaging over one day significantly improves the detection lim- 
its. This is consistent with the measured RV rms beeing domi- 
nated by the HF variations (see Table[T]). 

Finally, if we now consider set 1 data (no average), we get 
detection limits of 7, 8 and 16 Mj„p for periods of 5, 10, 100 
days, ie comparable to the limits found using the whole set of 
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Fig. 4. rms-based detection limits, for periods in the range 1- 
1000 days. Blue: all data, no average; red: all data, averaged over 
1 day; green: set 1 data, averaged over 1 day. 



data, when short periods are considered (limits for periods longer 
than typ. 200 days are not meaningful given the comparatively 
shorter duration of set 1). This similarity is due to the fact that 
the rms associated are comparable (see Table [T]). Finally, when 
averaging the data over one day, we find lower detection limits: 
3, 4 and 6.5 My„p for the same periods. These values are com- 
paratively better than those obtained after averaging the whole 
set of data, eventhough the time span and the number of data are 
smaller. This illustrates the usefulness of our adopted observing 
strategy, which naturally averages the stellar jitter. 

3.2. Periodogram-based detection limits 

Here we provide a brief description of alternative approaches 
we developped to take into account the structure of the stellar 
noise, and the associated results. We refer to a forthcoming pa- 
per (Meunier et al, 2012, in prep) for a detailed study of these ap- 
proaches, and of their respective merits, depending on the avail- 
able data. 

3.2.1. Highest peaks-based estimates 

For a given planetary mass and a given period, we compute the 
RV induced by the planet at the same observing dates as for the 
real data. We add this planet-induced RV signal to the observed 
one. We then identify the five highest peaks in the planet peri- 
odogram.0The mean power of the five highest peaks outside the 
window defined in the footnote is then compared to the power of 
the planet peak. For each planet mass and period, we simulate 
1000 realisations with different orbital phase between and In. 
The planet is declared detectable if in 99.9 percent of the simula- 
tions, the planet peak is higher than the mean of the five highest 
peaks. Given the steps adopted, the calculated detection limit has 
an accuracy of 0.5 Mj„p (set by the step taken in the computa- 
tionsfl 



- In practice, the simulated periodogram is affected by the injection 
of the planet not only at the planet period but also close to the planet 
period: several peaks, induced by the planet, are present close to the 
planet period. Therefore, we have to exclude these peaks; to do so, we 
compute the periodogram corresponding to a planet with period=Pp;„ 
for the same observing dates, and compute the associated median m and 
standard deviation s Meunier et al, 2012, in prep.). We then identify the 
peaks that are above a threshold of in+4s and located within 1 10% Pp/a 
from Ppi„ and we note Pf the period of the peak the furthest from Pp;„. 
The window size is then defined by min(10% Pp;„; Pf). 

^ Note: similar results were obtained using a criterion based on the fit 
of the RV data for periods corresponding to these five highest peaks. As 
the results are identical, we do not further describe this latter method in 
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Fig. 5. Detection limits found with the peak-based approach (5 
peaks, blue dashed/diamonds), and the correlation approach (red 
full line/stars). 



The detection limits obtained are shown in Figure|5] Typical 
values obtained for P ^ 5, 10, 100, and 500 days are resp. 2.4, 
3.8, 8.0, and 9.0 My„p (99.9% probability), hence in general sig- 
nificantly better than the rms-jitter-based ones when using un- 
averaged data; the difference is not so important when compared 
to the setl, rms-jitter-based data averaged over 1 day. This is 
consistent with the reduced jitter obtained once averaging setl 
data. For a 1000 day -period, the detection limit is larger than 50 

My,,/,. 

3.2.2. Correlation-based estimates 

In a second approach, we compare the periodogram of the planet 
alone to the periodogram corresponding to the total radial veloc- 
ities (ie the observed + planet induced RV). We compute the cor- 
relation between those two periodograms again for 1000 realisa- 
tions corresponding to the same planet mass, and the same planet 
period, and 1000 different orbital phase. In Figure|6] the obtained 
correlation values are shown as a function of the planet mass, 
for a given period (here 10 days). We see that for low planetary 
masses, the upper limits of the correlations for a given planet 
mass remain constant (ie do not depend on the planet mass). So 
do the maximum and minimum values at each mass. At these low 
masses, the planet is certainly not detected and so the correlation 
values there correspond to a 'no detection' -case. We define the 
detection limit as the lowest planet mass for which the small- 
est correlation value is higher than the threshold determined by 
the maximum correlation at the very low mass of 0.002 My„p. 
Hence the detection limit is the lowest mass for which 99.9% of 
the realizations have a correlation above that threshold. 

The detection limits obtained are shown in Figure|5] Typical 
values obtained are resp. 2.9, 1.8, 2.5, and 28.7 My„p for P ^ 
5, 10, 100, and 500 days (99.7% probability). Again they are 
better than the rms-based detection limits (no averaging), and 
more comparable (and actually slightly better) to those obtained 
by the highest peaks-based estimates. 




-0.5 



0.001 0.010 0.100 '.000 
Mass (Mjup) 



10.000 



the present paper but refer to Meunier et al, 2012, in prep, for such a 
description. 



Fig. 6. Correlations for 100 realisations at each planet mass ver- 
sus the planet mass (in My„p), for a period of 10 days. The solid 
line shows the mean for each mass. The dashed line indicates 
the threshold (determined from the maximum correlation at very 
low planet mass) and the arrow shows the detection limit. 



3.3. Planets with periods smaller than 1000 days: 
conclusions 

From the above study, we conclude that there is no planets with 
masses larger than 2.5 My„p and periods shorter than 100 days 
around /3 Pic . This excludes then the presence of massive hot 
Jupiters around yS Pic , whereas the star has certainly been sur- 
rounded by a massive protoplanetary disk and at least one outer 
planet has formed out of this disk. For planets with periods in 
the range 100-500 days, the detection limit is in the range 2.5- 
9 Mjiip. Finally, for planet-periods the 500-1000 day range, the 
present detection limit is in the brown dwarf domain. These are 
the best detection limits found around this star but also around 
any early type, rapidly rotating star, and these results clearly 
show the tremendous, yet sofar largely underestimated, possibil- 
ities of precise RV measurements of such stars. As far as /3 Pic 
is concerned, more RV data will help further constraining the 
detection limit at periods in the range 100-1000 days, and of 
course, also constraining the presence of planets at even larger 
periods. The present study shows that given the temporal struc- 
ture of the stellar noise, 2 hour-visits per year can in the future al- 
low to improve significantly the present detection limits for peri- 
ods > 200 days. Evenmore, planets with periods as large as 1000 
days or more could be tested. Coupled with other observing tech- 
nics such as interferometry jAbsil et al. 201 II I, aperture mask- 
ing (ILacour et al. 201 ll and AO imaging (Lagr ange et al. 20I0| l, 
this would allow to get for the first time a complete coverage of 
the mass-separation plane for giant planets and for separation as 
small as a few tenths of AU to several AU. 



3.4. Direct constrains on /3 Pic b mass 

The present data cover only 8 years, and do not therefore al- 
low to estimate detection limits for periods longer than typi- 
cally 1000 days, ie further than 2.5 AU. However, these data 
can be very usefull to constrain /3 Pic b mass as we have ad- 
ditional information on the planet position between 2003 and 
2011. According to present knowledge, the planet semi major 
axis is probably between 8 and 12 AU, with the most probable 
orbit at 9 AU. Its eccentricity is small (< 0.17) and its longi- 
tude of periastron (computed for instance from the line of sight) 
is not well constrained. Very importantly, in 2003, /3 Pic b was 
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close to a quadrature, moving away from us to pass behind the 
star, and in 20 11, close to the other quadrature, coming towards 
us ( ,Lagrange et al. 2010 ; Chauvin et al. 201 1 ). This is a very fa- 
vorable situation as we can expect (3 Pic RV variations induced 
by p Pic b to be close to the maximum amplitude at both dates. 
Knowing the companion positions in 2003 and 2011, using the 
fact that its orbital plane is seen edge-on, the RV data provide us 
with direct constrains on its mass, conversely to imaging which 
give model dependent (through brightness-mass relationships) 
estimates. 

To quantify an upper limit on [5 Pic b mass, we proceed as 
follows: 

- First, we measured the slope of the observed RV; the value 
found is 3.1e-5 km/s/day, with a fitting error of 1.4 e-6 
km/s/day. To estimate the actual error on the slope, we per- 
formed 1000 bootstraps of the data and measured the slopes 
obtained for each bootstrap. The average of the 1000 slopes 
is found to be 5.7 e-8 km/s/day (i.e. close to 0, as expected), 
and the dispersion is cr= 1.02e-5 km/s/day, a value unsurpris- 
ingly much larger than the error associated to the fitting pro- 
cess. We note that, as expected, this dispersion is identical to 
that obtained on 1000 realisations of random noise with an 
rms of 275m/s with the same calendar This shows that the 
measured slope is detected at a level of 3 sigma. 

- Then, we computed the detection limit for a planet with or- 
bital characteristics similar to those of yS Pic b. We assume 
that the planet was at quadrature beginning of January 2003, 
the currently most probable date dChauvin et al. 20lTb . We 
consider then a planet with a given semi-major axis (we as- 
sume a circular orbit), and a mass between 0. 1 and 25 M jup 
(steps of 0.1 My„p), and compute the RV induced by this 
planet on a 1.75 solar mass star, at the dates of the actual 
data. We add random noise with an rms of 275 m/s (corre- 
sponding to the jitter measured on [5 Pic RV data), and we 
perform a linear fit of these simulated RV, taking into ac- 
count the error bars measured at each date on /? Pic data, 
to estimate the slope of the variations together with associ- 
ated error bars. An example of the RV variations is showed 
in Figure |8] in the case of a 30 M y„p planet located on a 9 
AU (P=7400 d period): eventhough the sinusoidal shape of 
the variations is not obvious due to the high level of jitter, 
a significant positive drift is observed, which is not present 
in the p Pic RV (Figure [TJ. For each planet mass, we repeat 
this operation for 1000 different noise realisations in order 
to fully take into account the impact of the jitter noise. The 
values of the slopes thus obtained are showed in Figure|9]for 
a 9 AU orbit (in this figure, for the sake of clarity we have 
taken steps of 0.5 M y,,,, instead of Q.lMjup). We consider that 
a planet is detected if in 99% of the realisations, the slopes 
obtained are larger than 3xcr= 3xl.02e-5 km/s/dayQ 

- We then explored the impact of the date of quadrature, as- 
suming a range of -70;-i-130 days with respect to the most 
probable date, following lChauvin et al. 201T1 The impact on 
the planet mass is found to be 0.5-2 My„p, for semi-major 
axis in the range 8-12 AU. 

Table|2]give for semi-major axis of 8, 9, 10, 11, and 12 AU, the 
obtained upper limits. 

To try to improve these upper limits, we performed a simi- 
lar analysis using data averaged over 10 and 30 days. The upper 

^ Doing this, we assume that the level of dispersion in the slopes re- 
sulting from the boostrap of the yS Pic RV data is entirely due to stellar 
noise, which is a conservative estimation, as the observed slope indi- 
cates a drift at a 3 sigma level. 



Table 2. Constraints on p Pic b mass. See text. 



semi-major axis (AU) 8 9 10 11 12 
upper mass (Mj„,,) 9.6 12.0 15.4 20.0 24.8 



limits found are not improved which is due to the fact that be- 
cause of not optimal sampling of the data taken prior to 2008, 
the jitter is still high (100 m/s or more). 

In order to test the robustness of these results, as well as the 
impact of the planet eccentricity, we considered the statistical 
orbit distribution obtained from our MCMC analysis in Chauvin 
et al. (2011). For each model in this distribution, we computed 
the expected radial velocities at the observing dates, assuming 
a given planetary mass and eccentricities up to 0.17 following 
Cha uvin et al. 201T1 We then computed the statistical distribu- 
tion of the slopes after fitting the RV curve with a linear fit. The 
slope depends of course on the assumed mass of Beta Pic b. This 
allows us a compatible mass range, assuming that the measured 
slope must fall within the computed distribution; the slopes are 
then compared to the 3>cr limit. The results are found to be in very 
good agreement with those deduced from the previous analysis. 
This is illustrated in Fig [T] where the slope distribution is plotted 
for 3 masses (1.2, 9 and 24 MJup) and in the case of semi-major 
axis in the ranges 8.5-9.5 AU (left) and 1 1.5-12.5 AU (right). In 
the left figure, corresponding then to a ^ 9 AU, the distribution 
of slopes for a 24 MJup mass is well shifted compared to the 3cr 
limit, where as the one coiTesponding to a 9 MJup is closer. For 
a value of - 12 AU, the slope distribution with a 24 MJup mass 
is closer to the 3cr limit. 

Finally, we address the question whether these data also give 
constraints on the minimum mass of the planet. Assuming that 
the minimum slope is indeed equal to 3e-5 km/s/day, we can, 
as an exploratory exercice, derive a tentative lower limit to the 
mass of the planet if responsible for the slopes. Proceeding with 
a similar reasoning as above, and exploring again different dates 
of quadratures, we find lower masses in the range 1-2 MJup for 
semi-major axis between 8 and 12 AU. These conclusions are 
also compatible with the second approach described above, as 
distributions for e.g. 1.2 MJup are below the measured slope (see 
Fig. |7]i. The present values should be considered however only 
as an illustration of the potential of the present approach rather 
than definite results because the slope is so small and close to 
the limit; more data, with optimized temporal sampling are nec- 
essary to derive more reliable estimates of the p Pic minimum 
mass. 

3.5. Conclusion on p Pic b mass 

The present data have then allowed to constrain directly the p Pic 
b mass, to be less than 10, 12, 15.5, 20.0 and 25 MJup if orbiting 
respectively at 8, 9, 10, 1 1 and 12 AU. For the most probable or- 
bit, 9 AU, the upper limit is in the planetary mass regime, which 
means that definitively, the source detected is not a more massive 
object (brown dwarf). This is the first time an imaged extra-solar 
planet mass is directly constrained, and thus does not depend 
only on model-dependant brightness-mass relationships, p Pic b 
upper mass falls quite in the range predicted by the brightness- 
mass relations provided by the "hot start" models. On the con- 
trary, current cold-start models fail to explain a giant planet 
with the observed brightnesses at Ks ([Bonnefoy et al. 201 l| l, L' 
dLagran ge eTal. 20 10) and M (Quanz et al. 2010| ) bands. In the 
future, especially during the next coming years, new data with 
optimized sampling should allow to improve the mass range of 
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Fig. 7. Plots of the statistical distribution of the simulated radial velocity slope using the model developped in Chauvin et al (201 1), 
for three different masses : 1 .2 MJup (red line), 9 MJup (green line) and 24 MJup (blue line) and semi-major axis in the range 
8.5-8.5 AU (left) and 1 1.5-12.5 AU (right). The vertical grey/black bar shows the measured slope value, which is also the 3cr upper 
limit of the slope. 
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Fig. 8. RV induced by a 30 M j„p planet with a 7400 day-period 
(9 AU Sep.) and located at the position of /3 Pic b in 2003, when 
using the actual data sampling: grey crosses: without jitter noise, 
black squares: with a noise of 275 m/s. The grey line indicates 
the slope derived from the linear fitting of RV variations, and the 
black line indicates the slope of /3 Pic real data deduced from 
linear fitting. 



jS Pic b. The present combination of RV and high contrast data on 
exoplanets is a pioniering one; it will be more generalized in the 
future, when more young planets will be imaged at short separa- 
tions from their stars, and will give thus a precious opportunity 
to test brightness-mass predictions. 
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